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Taking up and extending earlier suggestions, we show how two- and threedimensional shapes of second- 
order HBT correlations can be described in a multivariate Edgeworth expansion around gaussian ellipsoids, with 
expansion coefficients, identified as the cumulants of pair momentum difference q, acting as shape parameters. 
Off-diagonal terms dominate both the character and magnitude of shapes. Cumulants can be measured directly 
and so the shape analysis has no need for fitting. 
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I. INTRODUCTION 

Early measurements of the Hanbury-Brown Twiss (HBT) 
effect made use of momentum differences in one dimension, 
for example the four- momentum difference ^inv [l]- The 
huge experimental statistics now available permits measure- 
ment of the effect in the three-dimensional space of vector 
momentum differences q = pi — P2 and, in many cases, also 
its dependence on the average pair momentum K = (pi + 
P2) /2. Increasing attention has therefore been paid in the last 
decade to the second-order correlation function in its full six- 
dimensional form. 



more traditional continuous version 111211 with a gaussian ref- 
erence /o- For the purpose of analysing the shape of the ex- 
perimental coiTelation function in HBT, we hence define the 
measured nongaussian probability density as 



/(q,K) 



^2(q,K) 
JdqR2{q,K)- 



(3) 



where ^2(q) is itself a normalised cumulant of pair counts 
ifisll . For the reference distribution (null case), we take the 
multivariate gaussian, which in its most general form is 



C2(q,K) - 1 + R2iq,K) 



(1) 



/o(q,K) 



exp [- 1 {qi - A-,-) (A- ' ),7 (qj - Xj)] 
(27i)0/2(detA.)i/2 



(4) 



J dpi dp2 p(pi,P2)5(pi -p2-q)5(^(pi +p2)-K) 
J dpi dp2 P'''=^(pi ,P2) 5(pi - P2 - q) 5(^(pi + P2) - K) ' 

with p the density of like-sign pairs in sibling events and 
p"^^ the event-mixing reference. While C2 data can be visu- 
alised and quantified reasonably well in two dimensions |l2|,|3t], 
it is harder to quantify three- or higher-dimensional correla- 
tions. Projections onto marginal distributions are inadequate 
0\, while sets of conditional distributions ("slices") require 
many plots and miss cross-slice features. 

Under these circumstances, efforts to quantify the shape of 
the multidimensional correlation function with Edgeworth ex- 
pansions |5, 6] or spherical harmonics |4] represent welcome 
progress. More ambitious programmes seek to extend con- 
nections between gaussian source functions and the "radius 
parameters" of the correlation function to sets of higher-order 
coefficients using imaging techniques 10. H. 1^ and cartesian 
harmonics Llftl . 

In this contribution, we extend the Edgeworth expansion 
solution proposed in (?,^ to a fully multivariate form, in- 
cluding cross terms. Generically, the intention is to expand 
a measured normalised probability density /(q) in terms of a 
reference density /o (q) and its derivatives, 

/(q) = /o(q) { Edgeworth expansion in q } , (2) 

so as to characterise /(q) by its expansion coefficients. While 
we have previously made use of a discrete multivariate Edge- 
worth form with poissonian reference /o to describe multi- 
plicity distributions ifTlll . the shape analysis of R2 requires the 



where D is the dimensionality of q, Einstein summation con- 
vention is used (here and throughout this paper), the are the 
first "q-cumulants" of /o ifliilTsll . 



?i,(K)=/flfq/o(q,K)^,-, 



(5) 



A-,j(K) is the covariance matrix (the set of second-order q- 
cumulants) in the components of q, 

= {Idqfoqiqj)~{Jdqfoqi){fdqfoqj), (6) 

and Xjj^ the inverse matrix. While we suppress K in our no- 
tation from now on, all results are valid for K-dependent first 
moments A,, and covariance matrix elements Xij. 

II. REFERENCE DISTRIBUTION 

The vector difference is normally decomposed into com- 
ponents q = {qi,q2,q3) — {qo,qs,qi) in the usual (out, side, 
long) coordinate system; for illustrative purposes we will also 
make use of a two-dimensional vector q = (^1,^2)- (This 
is not the two-dimensional decomposition into {qt,qi) used 
in some experimental HBT analyses 13, Hd, because qt = 
ill + I's)^^^ is always positive, while (^1,^2) in our two- 
dimensional example can be positive or negative.) For the 
two-dimensional decomposition, the covariance matrix 



X = 



All ^12 
A12 A22 



(7) 
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has the inverse 



mensions the covariance matrix 




(8) 



where we have introduced standard deviations a,- 



well as the Pearson correlation coefficient p 
^i2/VA.iiA,22 andx= 1 - Iil2>,l6.1 



12/(0102) = 
Similarly, in three di- 



1 = 



^ K>o K>s K)i \ 

\).<i ^SS \<il 

\ Ka Ki hi J 



(9) 



has the general inverse 



detA, 



i o2o2(l-p 



detA, 



-QoOsO] (Pf«-p„/P,/) 
\ -0„0?a; (p„/ - p,„Ps/) -O^OjO/ (p.v/ - p,«p„/) 



-0„a,0^ (Po,. - P„/P,/) -OoO^O/ (po/ - Po,p,/) ^ 

-o^o.,a/ (Pi/-Pf„Po/) 



where p,^ = djj {(5i<3j) and the determinant is given by 



2 2 2 222/222 \ 

detA = Xoohshl—^o'^sl^h'is'^ol^hlhs + '^h'sKil^sl = 0„OjO/ (1 — Pj/ — Po/ — pf„ + 2p„sP„/Ps/) 

I 



(10) 



(11) 



(12) 



For azimuthally symmetric sources [17], po 
the inverse simplifies to 



A-' = 



1 



-Pol 











-Pol 



1 



/ 2Ri, 2Rl, \ 



- Psi = 0, so that fji — J dqf{q) qt, and — J dqf{q) qigj, and in general 

/Jijk... = dqf{q)qiqjqk... . (14) 



2Rl 






24 J 



(13) 



Identifying in the second part of Eq. (l3[ the inverse cu- 
mulant matrix with the usual radii Rfj of the parametrisation 

/o ~ exp[— we note that the notation is mis- 
leading in that a positive covariance between the out and long 
directions, p^/ > 0, results in a negative R^^j. 



III. MULTIVARIATE EDGEWORTH EXPANSION 
A. Derivation 

In order to derive the Edgeworth expansion, we need to dis- 
tinguish between the moments and cumulants of /o(q) and 
/(q) respectively. The cumulants of the reference /o have 
been fully specified already: the order- 1 and 2 cumulants are 
the set of (initially free) parameters A, and A,j respectively, 
while all cumulants of order 3 or higher vanish identically 
lfl2ll for the gaussian reference (01). For the measured nongaus- 
sian /(q), we denote the first- and second-order moments as 



Cumulants K,jit... of /(q) are found from these moments by 
inverting the generic moment-cumulant relations lfl2ll 



Kij + K,K j , 

Kijk + KiKjk + KjKki + KkKij 



(15) 
(16) 

(17) 



where we have introduced the notation [3] to indicate the num- 
ber of index partitions, and therefore terms, of a given combi- 
nation of k's. The relations of order 4, 5 and 6, which we will 
need in a moment, are 

t^ijkl Kijkl+KiKjkl[4]+KijKkl[3] (18) 
+ K,K;Ki./[6] +K,-K;KiK/, 
l^ijklm = Kijklm+KiKjklm[5]+KijKki„,[lO] (19) 

+ K,K;K^K,„ [10] + K,KjKi,K/K„, , 

l^ijklmn = l^iylZmn + l^i'l^ jWmn [6] + K,jK/;/„,„ [15] (20) 

+ KiKjKklmn [15] + K,yi;K,,„„ [ 1 0] , +K,KjiK,,„„ [60] 

+ K,K^Ki,K/„„, [20] + K;,-Ki,/K,„„ [15] + K,K;Ki.,K™ [45] 

+ K,K^Ki,K/K,„„ [15] + K,K^Ki.K/K,„K„ . 
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For identical particles, all moments and cumulants are fully 
symmetric under index permutation. 

The derivation of the Edgeworth expansion starts with the 
generic Gram-Charlier series lfl2[[l8l[l9ll . which is expressed 
in terms of differences between the measured and reference 
cumulants 



r|, = K,-A.,, 
'nijk = ^ijk-hjk etc., 



(21) 
(22) 
(23) 



and moment-like entities t^tjk... which are related to the T\ijk... 

i.e. 



by the same moment-cumulant relations (fT5]l-(l20li. 

C; = r|; = K; - A.;, 

C<7 - llu + m^j = (k,; - hj) + (K; - h) (K; - 
and so on. The Gram-Charlier series 
/(q) 



(24) 
(25) 



/o(q) 



(26) 



is an expansion in terms of the ^s and partial derivatives 

1 aVo 



hijiq) 



hijkiq) 



/o dqidqj ' 

1 aVo 



(27) 



(28) 



/o dqidqjdqk 



etc., 



(29) 



which for gaussian /o are called hermite tensors; they will be 
discussed below. 

The generic Gram-Charlier series is reduced to a simpler 
HBT Edgeworth series in three steps. First, the freedom of 
choice for the parameters A,, and A,,^ of the reference distri- 
bution (|4|i allows us to set these to the values obtained from 
the measured distribution; i.e. we are free to set A,- = K, and 
Xij = Kij, SO that r\i = r\ij = and hence ^, = = 0. 

Second, we make use of the fact that all cumulants of order 
3 or higher are identically zero for the gaussian distribution, 

A,7A- = XijU = • • • = 0, so that Z^ijk = Kijk, ^,i}kl = y^ijkl, ^ijklm = 

Kjjkin, and in sixth order !^ijki,„„ = Yiijkimn + ^ijk^imn [10] . 

Finally, the contribution to the correlation function C2(q) of 
the momentum difference q^p = Pa=i — Pp=2 of a given pair 
of identical particles (a, P) is always balanced by an identical 
but opposite contribution qp„ = pp^i — pa=2 = ^ qap by the 
same pair, so that C2(q) must be exactly symmetric under "q- 
parity". 



C2(-q)=C2(q). 



(30) 



This implies that all moments and cumulants of odd order of 
the measured /(q) must be identically zero, Kijk = "i^ijUm = 0, 
so that terms of third and fifth order and the KijkKi^n contribu- 
tion to sixth order are also eliminated. 

The end result of these simplifications is a multivariate 
Edgeworth series in which only terms of fourth and sixth order 
survive. 



/(q) 
/o(q) 



= 1 + ^ K^ijkl hijkl (q) + 5! ^ijklmn hjjklmn (q) - 



(31) 



For three-dimensional q, there are 8 1 terms in the fourth-order 
sum and 729 in sixth order, but due to the symmetry of both 
the K and h, many of these are the same. Defining n — n\ + 
ni + "3, we introduce the "occupation number" notation 



(-iravo 



1 



/o (3^i)"i(a<72)"2(3^3)"3' 



(32) 



and correspondingly define cumulants Cnyn2n^ as K,|/2...,,j with 
n\ occurrences of the index 1, n2 occurrences of 2, and n3 oc- 
currences of 3 in {i\ ■ --in), e.g. Cm = K1223 = K3122 = . . .. 
Similar definitions hold for //„,„, and C„,„2 for the two- 
dimensional case. Combining terms in (IJTT i. we obtain for 
the two- and three-dimensional cases respectively. 



/(q) 
/o(q) 



= 1 + ^ {C40//40P] +4C31//31 [2] +6C22H22} 

+ ^{C60//60[2]+6C5li/5l[2] 



/(q) 

/o(q) 



+ 15C42//42[2]+20C33//33 } + ... (33) 

1 + ir {Qoo^40o[3] +4C31o//31o[6] 

+ 6C220i^220[3] + 12C21l//21l[3] } 

+ ^ {Ceoo^/^soop] +6C5ioH5io[6] 

+ 15C42oi?42o[6] +3OC411//411 [3] 
+ 20C330//330[3] + 6OC321//321 [6] 
+ 9OC222//222} + ... (34) 

where the square brackets here indicate the number of dis- 
tinct cumulants related by index permutation to those shown. 
In two dimensions, we therefore have 5 distinct cumulants of 
fourth order and 7 of sixth order, while in three dimensions 
there are 15 distinct fourth-order and 28 sixth-order cumulants 
respectively. We note that these cumulants can be nonzero 
even when the reference gaussian is uncorrected, i.e. even if 
the Pearson coefficients are zero. 

B. Hermite tensors 

In the Edgeworth series ( |33] ) and (l34l i. the cumulants C are 
coefficients fixed by direct measurement, while the hermite 
tensors H are, through eqs. (l27])-(|29b and explicit derivatives 
of (|4|i, known functions of q. Defining dimensionless vari- 
ables 



Zi = 



qi_ 



(35) 



which can also be written in terms of the usual radii as 
Zi — V^qiRii, the lowest-order hermite tensors are, for the az- 
imuthally symmetric out-side-long system. 



hi = HiQo 
hi =Hoio 
hj, — Hqqi 



Zl - pZ3 
Z2_ 

Z3 - PZI 

5C03 



(36) 
(37) 
(38) 
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FIG. 1: Surface plots of Gaussians times individual hermite tensors, 
fnHnim in two dimensions. Tensors //o4 and //13 are images of //40 
and Hi,\ mirrored through the zi = zi diagonal. Axis labels are in 
units of \/2o/. 



Fourth-order derivatives ( [32] l of the gaussian (|4|i yield, 



Haqo = /it-6/i?X-/+3V/Vi', (39) 

Horn = /i3 — 6/13X33' + 3X33^X33', (41) 

//310 = h\h2-3hxh2X^l, (42) 

i/l30 = /!^/!I-3/!l/l2A.22\ (43) 

H013 = /ii/z2- 3/12/13X3-3', (44) 

i/03l = /^^/^3-3/^2/^3X22^ (45) 

i/301 = /z^/i3-3/zi/!3Vi'-3/!?A.i3' + 3A.7i'A.i3', (46) 

//i03 = /i3/ii-3/!i/i3A.33'- 3/1^X13' + 3A.33'A.j3i, (47) 

H220 = h\hl^h\l22 -hlKi +Kl^22^ (48) 

^^022 = h^h2 — h^}^22 ^^2^33 ^ '^22 '^33 ^ (49) 

H202 = h]hj-h]X^^ -hjX];! -4hih3X-^ 

+ A,j j'A-33' + 2A-j3'A,j3' , (50) 



//211 = h\h2h3-2h\h2'k\^ -h2h3X\^ , (51) 

^112 = /ll/l2/l3 — 2/Z2/J3A-J3' — /zi/Z2A-33', (52) 

//i2i = hihlh3-hlX-^^ -hih3X22 +K2K3: (53) 

where the inverse cumulant elements Xjj ' are functions of the 
parameters A-; and Xij that can be read off from Eq. ( fT3] ). The 
differences between various permutations of (nin2«3) above 
arise from the fact that A-jj' = A-23' = for azimuthal symme- 
try. 

Note that only one of the above hermite tensors can be writ- 
ten in terms of hermite polynomials at this level of generality, 
namely //o4o = H4 {z2 )/02- Generally, the hermite tensors fac- 
torise into products of hermite polynomials H„ (z, ) only if all 
Pearson coefficients p,^- in the gaussian reference are zero, 

Hn.n.n, (p,7=0) = fl (54) 

7=1 

In sixth order, the tensors are generically 

hijkimn ~ hih jhkhihmhn — hih jhkhiX^l [15] 

+ h,hjX^X.n [45] - ^ijXX'n [15] (55) 

where again the square brackets indicate the number of dis- 
tinct index partitions. Sixth-order tensors //„,„2«3 can then be 
constructed from these as usual, for example 

HeoQ = h^l-lShlX-l+ASh^lf-lSiXilf (56) 

which closely resembles the Hermite polynomial //g (z) = z*' — 
15z^ + 45z^ — 15 but reduces to the latter only when pd — 
and hence X^^ = 

C. A gallery of shapes 

In Fig. 1, we show surface plots for individual fourth-order 
hermite tensors times the two-dimensional reference gaussian, 
foH„^„^, with p set to zero. As these are plotted in terms of 
V2zi — \/2qi/<5i lEoll . the axes are scaled by the standard de- 
viations, meaning that all gaussians with p = will be circular 
in (zi ,22) plots. The individual hermite tensors clearly reflect 
the symmetry of their respective occupation number indices n,- 
and probe different parts of the (zi ,Z2) phase space as shown. 
Comparing the fourth-order terms of Fig. 1 with the sixth- 
order ones of Fig. 2, we note that the latter probe regions up 
to several V2c?;- 

In order to exhibit the influence of combinatoric factors, we 
show in Fig. 3 individual terms of the two-dimensional Edge- 
worth series in the form /o(q)(l +/V!i«2C„i«2^«i«2(<j))' 
where the combinatoric factors F,^„^ are fixed in (l33T l. All 
fourth- and sixth-order cumulants C„,„2 have been set to 1.0 
and 2.0 respectively. (This is obviously for illustrative pur- 
poses only; in real data, smaller values are expected and 
shapes will be more gaussian than those shown here.) The 
plots for H31 and H\3 illustrate the correspondence between 
index permutation and symmetry about the zi = Z2 axis. We 
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FIG. 2: Contour plots of Gaussians times individual sixth-order her- 
mite tensors, foHi^,,-, in two dimensions. Red-green (light grey) ar- 
eas represent hills while blue-red (dark grey) areas are valleys. Note 
how regions of phase space at a distance of several o, from the peak 
are probed. 

note that the diagonal terms H„q have little influence on the 
overall shape, while the off-diagional ones have a larger ef- 
fect, not least because of the combinatoric prefactors. 

Testing the influence of fourth- versus sixth-order terms, we 
show in Fig. 4 some "partial" two-dimensional Edgeworth se- 
ries including only fourth-order terms, only sixth-order terms, 
and both orders; again, cumulants were set to the arbitrary 
values of 1 and 2 respectively. 

In Fig. 5, a selection of shapes for individual terms 

/o(q)(l + F„|„2„3C„,„2»3^/^p,i«2n3) the three-dimensional 
Edgeworth expansion ( [34l l is shown. While in the two- 
dimensional case full contour plots could be shown, the sur- 
faces shown here in each case represent only a single con- 
tour. In Fig. 6, we show two examples with two selections of 
fourth-order cumulants nonzero; the shape obviously depends 
strongly on their selection and magnitude. Clearly, effects of 
the different cumulants on the overall shape often cancel out. 
We emphasize again that the shapes shown are for illustrative 
purposes only and do not represent real data. 




H42 «33 



FIG. 3: Contour plots of some individual terms in the Edgeworth ex- 
pansion ([33} of the form /o(l + F„^„^C„^n2Hnin2) Arbitrary and unre- 
alistically large values of C„,„, = 1 for fourth order and C„,„, = 2 for 
sixth order were chosen, while the combinatoric prefactors F„ ■ 11-, are 
fixed in ( I33t . Both H^i and //13 are shown to exhibit their symmetry 
about the zi = Z2 diagonal. 



IV. DISCUSSION 

The multivariate Edgeworth expansions (l33Tl-(l34li appear 
to be a promising tool for quantitative shape analysis in HBT. 
While the real test will be to gauge their performance in ac- 
tual data analysis, they do seem to have the right features and 
behaviour A number of issues deserve further comment: 



1. It has been noted previously ifTilfTsIl that the traditional 
radii of a gaussian-shaped ^2(q) could be found by di- 
rect measurement rather than from fits. In the present 
formulation, this amounts to the direct measurement of 
the second-order cumulants A,,,-, which can be directly 
converted to "radius parameter" form via Eq. ( fTSl ). Go- 
ing beyond Refs. iJJlllJ], we suggest that higher-order 
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(c) (f) 

FIG. 4: Combined terms for the two-dimensional Edgeworth: fourthi- 
order terms only in (a) and (d); sixth-order terms only in (b) and 
(e); and all terms in (c) and (f). Fourth-order cumulants were set 
arbitrarily to 1 and sixth-order ones to 2. In (a)-(c), all cumulants 
are equal and nonzero, while in (d)-(f), cumulants C04, C13, C06, 
Ci5 and C24 were set to zero in order to illustrate the possibility of 
asymmetric shapes. 



cumulants can be measured directly also. 

2. Many people have rightly expressed concern that these 
radii do not adequately represent the true shapes and be- 
haviour of HBT correlations. Our Edgeworth expansion 
confirms that such radii are clearly not the whole story, 
but that they do represent the appropriate lowest-order 
approximation (for gaussian reference) with respect to 
which nongaussian shapes should be measured. 

3. We have demonstrated that it is imperative to write 
Edgeworth expansions in a fully multivariate way: the 
combinatoric prefactors F„^„^„^ in ( l34l i are large for 
multivariate "off-diagonal" cumulants, while the influ- 
ence of diagonal cumulants is strongly suppressed due 
to their small prefactors. The cumulant C211, for exam- 
ple, has a weight 12 times larger than C400, and indeed 




FIG. 5: Typical contours for three-dimensional terms /o(l + 
n2ni ^«in2«3 ) with indices as shown. Note that only a single 
contour is shown in each case. 




FIG. 6: (a) Typical contour plot for combined Edgeworth with cu- 
mulants C310, C220> C211 and their permutations set to 0.3 with all 
others (including 6th order) set to zero, (b) Combined Edgeworth 
with cumulants C22Q, Qoi^ Co02> C21I' C112 set to 0.15 and all 
others to zero. Only single contours are shown. 



the entire expansion is dominated by the off-diagonal 
cumulants. Furthermore, even large diagonal cumulants 
do not change the shape much, as a glance at Fig. 3 will 
confirm. 

4. Deviations from gaussian shapes are consistently quan- 
tified by the sign and magnitude of higher-order cumu- 
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lants, which are identically zero for a null-case pure 
gaussian /(q). The Edgeworth expansion using these 
cumulants, while recreating the shape of /(q), there- 
fore at the same time provides a quantitative framework 
for comparison of different shapes. 

5. Operationally, we suggest a procedure of successive ap- 
proximation, whereby in a first step all elements of the 
covariance matrix K,j — Xjj are measured, thereby de- 
termining all the a's and the Pearson coefficient; this 
is equivalent to the usual determination of radii. This 
is followed by measurement of the set of fourth-order 
cumulants C,^,n„^. The measured numbers for fourth- 
order cumulants then represent the basis for shape quan- 
tification and comparison. If statistics permit, sixth- 
order cumulants can be added as a further refinement. 

6. The q-cumulants proposed here are numbers rather than 
functions of q. From the viewpoint of compactness of 
description, this will be an advantage compared to the 
shape decompositions in terms of spherical and carte- 
sian harmonics ESl, in which each coefficient is a 
function of |q|. It may, however, in some cases be bet- 
ter to see the detail provided by such functions. 

7. The procedure outlined above involves no fits. This rep- 
resents a major advantage over fit-based quantification 
in two ways: 

Firstly: In three-dimensional analysis, typical fits are 
dominated by phase space, i.e. by the fact that there are 
many more bins at intermediate and large |q| than at 
small |q|. This dominance suppresses the influence of 
the most interesting region on the for best-fit values 
of the parameters. In Ref. [311, for example, we found 
that the regions of intermediate |q| dominated the shape 
and quality of various fits. 

Secondly, as shown in Fig. 2, cumulants are sensitive to 
the tails of distributions, and they will hence access the 
same information as these fits and parametrisations, but 
in a more direct and sensitive way. It is well known that 
a direct fit to a probability distribution that is close to 
gaussian is an ineffective and inaccurate way to quan- 
tify nongaussian deviations, while cumulants do so in 
the most direct way possible. 

8. It remains to be seen how the proposed procedure fares 
when the practical experimental difficulties of finding 
/(q) and the higher-order cumulants come into play. 
Much will also depend on the size and accuracy of 
statistical errors. Fortunately, current sample sizes are 
large enough to warrant some optimism in this respect. 

9. The traditional chaoticity parameter X remains unde- 
termined within the present Edgeworth framework, be- 
cause it cancels already in the definition Q of /(q). For 
a given level of approximation (gaussian only, fourth- 
order cumulants, sixth-order), it and the overall normal- 
isation factor y may be recovered afterwards by using 



(l34b in a two-parameter fit mode using parametrisation 

C2(q) = y [1 + A-/o(q) (Edgeworth expansion) ] 

with the previously experimentally-determined radii 
and C„j„2«3 treated as constants, with y and X the fit pa- 
rameters. 

10. We note the importance of the parity argument 
C2(— q) = C2(q) in eliminating odd-order terms in the 
Edgeworth expansion. The parity argument falls away, 
however, in variables where this symmetry does not 
arise; for example, any one-dimensional Edgeworth ex- 
pansion involving only positive differences (e.g. in qi„y) 
would have to include third- and fifth-order terms. 

A corollary of the parity argument is that three- 
dimensional correlations may not be represented in 
terms of positive absolute values of the components 
{qo,qs,qi) as this destroys the underlying symmetries. 
The best one can do to improve statistics is to combine 
bins that map onto each other under the transformation 
q ^ — q and thereby eliminate four of the eight octants 
in the three-dimensional (^o,^s,^/)-space. 

11 . In the present formulation, any dependence on average 
pair momentum K resides in the cumulants: all K/yj^ . , 
including the traditional radii and the Pearson coeffi- 
cient, are in principle functions of K. 

12. The Edgeworth analysis set out in this contribution is 
based on a gaussian reference /q. Shapes that differ 
significantly from gaussian will not be described well 
in either the Edgeworth framework or the spherical or 
cartesian harmonics frameworks. One should not, for 
example, expect power laws such as a pure Coulomb 
wavefunction (whose square tails off like |q|^^) to work 
in a gaussian-based Edgeworth expansion. Indeed, it 
is known that large cumulants can lead to a situation 
where the truncated Edgeworth expansion of /(q) be- 
comes negative in some regions. It is therefore suitable 
only for shapes that do not deviate strongly from gaus- 
sians; for strong deviations, other expansions will be- 
come a necessity. 

13. The Edgeworth framework is easily extended to the 
case of nonidentical particles. In that case, cumulants 
of all orders will have to be measured. It may well 
be that the fluctuations of lower-order quantities render 
the measurement of higher-order cumulants impossible, 
and great care will clearly have to be taken. 
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